Unifying several separability conditions using the covariance matrix criterion 
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We present a framework for deciding whether a quantum state is separable or entangled using 
covariance matrices of locally measurable observables. This leads to the covariance matrix criterion 
as a general separability criterion. We demonstrate that this criterion allows to detect many states 
where the familiar criterion of the positivity of the partial transpose fails. It turns out that a large 
number of criteria which have been proposed to complement the positive partial transpose criterion 
- the computable cross norm or realignment criterion, the criterion based on local uncertainty 
relations, criteria derived from extensions of the realignment map, and others - are in fact a corollary 
of the covariance matrix criterion. 

PACS numbers: 03.67.-a, 03.65.Ud 



I. INTRODUCTION 

Entanglement is the feature of quantum theory that 
renders it crucially different from a classical statistical 
theory. It also plays the central role in quantum infor- 
mation science, as a resource for information processing 
tasks. Consequently, a lot of effort has been made in the 
last decade to understand the meaning and the structure 
of entangled states [l], ■ One of the most elementary yet 
notorious questions is how to find good criteria to decide 
whether a state is entangled or classically correlated in 
the first place. Formally, one asks whether a given state 
p in a bipartite system is contained in the convex hull of 
product states and can hence be written as 
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where the coefficients pk form a probability distribution. 
If so, all correlations can come from classical shared ran- 
domness, and a state is called classically correlated or 
separable. Otherwise, p is an entangled state. The deci- 
sion problem of deciding whether a state is entangled or 
separable is known to be a computationally hard prob- 
lem in the physical dimension, for a certain scaling of the 
error in the weak membership problem Yet, the prob- 
lem one typically faces is the one where one has just a 
physical state given - having some fixed dimension - and 
one aims at finding criteria to make a judgment based 
on these criteria. Not surprisingly, given the central sta- 
tus of entanglement in quantum information theory, a 
lot of effort has been devoted to identifying such good 
and practical and computable criteria for separability in 
composite quantum states. 

Historically, the first criterion of this sort was derived 
from the observation that every separable state will have 
a positive partial transpose, and that the positivity of 
the latter can hence be used as an entanglement crite- 



rion (PPT criterion) 0, Q. This criterion later turned 
out to be necessary and sufficient for separability for low 
dimensional systems (2x2 and 2 x 3), whereas in higher 
dimensions this is no longer the case [||. The PPT crite- 
rion is an example of a criterion based on positive maps: 
In fact, it has been proven that a state p is separable if 
and only if for any positive map A the operator (l(g>A)(p) 
is also positive [5j. Consequently, the systematic inves- 
tigation of positive maps has led to a number of new 
separability criteria Q. A quite remarkable criterion of 
this type is the reduction criterion, which is equivalent to 
the PPT criterion for 2 x 2 and 2x3 cases and weaker for 
higher dimensions Q. There are also other criteria that 
turned out to be directly related to the PPT criterion: 
The majorization criterion Q and also entropic criteria 
10] have been shown to be weaker than the PPT criterion 
11,0]. Moreover, one can extend the PPT condition to 
a test based on a complete hierarchy of symmetric exten- 
sions, where each step constitutes a semidefinite program 
[l3l ] (another complete family of semidefinite tests has 
been described in [T3]). In such a hierarchy, every en- 
tangled state is necessarily detected as such in some step 
of the hierarchy. However, the steps in the hierarchy re- 
quire more and more computational effort, which makes 
this approach difficult already for modest system sizes. 

Apart from these criteria that are directly related to 
the PPT criterion, a number of other separability crite- 
ria have been suggested where such a connection seem- 
ingly does not exist. The most prominent criterion of 
this type is the computable cross-norm or realignment 
criterion (CCNR) [151 ] . Other criteria in this cate gory 
use the Bloch representation of density matrices (lH Il7j. 
local uncertainty relations (LURs) [181 ]. local orthogonal 
observables (l9| . or extensions of the realignment map 
[13, Hl|- These criteria are complementing the PPT cri- 
terion in an interesting way: In fact they detect some 
states as being entangled where the PPT criterion fails. 
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They also do not rely on positive maps. 

At first sight, one might think that these criteria form 
a collection of quite beautiful, but strangely disconnected 
results. They have been derived using a variety of unre- 
lated methods, and their connection often seems quite 
unclear. It is the main purpose of this work to develop a 
framework for the systematic understanding of all these 
latter approaches. 

In Ref. [13] we have proposed to investigate the separa- 
bility problem using covariance matrices (CMs) of certain 
observables. In this context we have developed a sepa- 
rability criterion in terms of covariance matrices (covari- 
ance matrix criterion or CMC). We have shown that the 
CMC is, when augmented with appropriate local filter- 
ing, despite its simplicity a surprisingly strong entangle- 
ment criterion, which can detect states where the PPT 
criterion fails and which is at the same time necessary and 
sufficient for two qubits. Here we complete this approach 
and present new results in various directions. Specifically, 
we show that a number of criteria which have been pro- 
posed to improve the PPT criterion - namely the CCNR 
criterion [15j . a criterion using the Bloch representation 
16, 03, the LURS [HI], and recent criteria from Refs. 
2fJ, |21| - follow directly from the CMC. In this way the 
CMC can be seen as complementary to the PPT crite- 
rion. We also tighten previous formulations of the CMC. 
We discuss several examples, and compare the perfor- 
mance of the criteria to instances of random states from 
a families of bound entangled states. 

This manuscript is organized as follows: In the second 
section, we introduce CMs and discuss their mathemati- 
cal properties. Those readers who are mainly interested 
in separability criteria may only consume Definitions III. II 
and III.2I and Propositions III. 71 III. 81 and III. 121 and may 
then directly jump to Section III. In that Section, we 
introduce the CMC and evaluate it in several different 
ways. By doing this we establish the mentioned connec- 
tion to the other separability criteria which will turn out 
to be corollaries of the former. In Section IV of the pa- 
per we will consider the connection between the CMC 
and the LURs. In the fifth section we will scrutinize the 
CMC for the two qubit case. In Section VI we will assess 
the strength of the mentioned criteria by considering a 
family of bound entangled states. We will then conclude 
and elaborate on possible extensions of the work pre- 
sented here. Some more technical proofs of our theorems 
will finally be presented in the Appendix. 



II. COVARIANCE MATRICES 

In this section we will investigate covariance matrices 
as our main tool. In the first subsection we will introduce 
the different definitions of CMs [U [H [H| and fix our 
notation. In the second subsection we will address the 
question to which extent CMs can be used as a unique 
description of quantum states besides density matrices. 
Finally, in the third and fourth subsection we will men- 



tion and prove some useful properties of CMs, which will 
be used later in our study of entanglement. 



A. Definition of covariance matrices 

In what follows let p be a pure or mixed quantum 
state, described by a (positive) density operator in a d- 
dimensional Hilbcrt space TC and let {Mk : k = 1, . . . , N} 
a suitable set of observables. Unless stated otherwise, we 
will always assume that these observables are orthonor- 
mal observables with respect to the Hilbert-Schmidt 
scalar product between observables, i.e., they fulfill 



tx(MiMj 



(2) 



Furthermore, we will typically assume that the Mi form 
a complete basis and span the whole observable algebra. 
This implies that there are N = d 2 different Mi, and 
that any other observable can be expressed as a linear 
combination of the Mi. 

As an example for such a set of observables for the 
case of a single qubit, one can consider the (appropriately 
normalized) Pauli matrices, 



Mi 



1 



M 2 



M 3 



y 



Ma 



V2' 



(3) 



12' ' V2' V2 
We can now formulate the main definitions for this work 



Definition II. 1 (Covariance matrix). The d 2 x d 2 co- 
variance matrix 7 = j(p, {Mk}) and the d 2 x d 2 sym- 
metrized covariance matrix 7 s = 7 (p, {Mk}) are defined 
by their matrix entries as 



'hi 



{MiMj) — (Mi)(Mj) 
(MiMj) + (MjMi) 



{Mi)(Mj 



(4) 
(5) 



Sometimes, the difference between the linear part of a CM 
and the nonlinear part becomes relevant. Therefore, we 
define the linear part 0/7 as Qi j = (MiMj) and the linear 
part of the symmetric CM as gfj = (MiM 3 + MjMi)/2. 

We will often for simplicity of notation also write "f(p) 
or 7({Mfc}) instead of j(p, {Mk}), or simply 7. We will 
also sometimes indicate with respect to what state an ex- 
pectation value is taken, so (Mi) = (Mi) p . It is straight- 
forward to see that 7 is a complex Hermitian matrix. The 
matrix 7 s in turn is real and symmetric. Both 7 and 7 s 
are positive semidefinite, 7,7 s > (24j. 

Note finally that for odd d, there is another basis of 
orthonormal observables that can equally be used and 
that is commonly employed in the mathematical physics 
literature in the context of discrete Weyl systems [25| . 
Let A(0, 0) be the parity operator that maps P(0, 0) : 
\x) 1 — > |— x), where \x) <G {|0), . . . , \d— 1)}, meant modulo 
d. Then, for (q,p) G 7L 2 d let 



P(q,p) =W(q,p)P(0,0)W(q, P y 



(6) 
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the translated versions of P(0, 0) in discrete phase space, 
where W(q,p) are the discrete Weyl operators [26]. The 
operators {Mr qjP \} = {P(q,p)Vd} then form a set of 
Hilbcrt-Schmidt orthonormal Hermitian matrices. This 
is the standard set of observables when phase-space 
methods are made use of. 



B. Covariance matrices for bipartite systems 

In the focus of this work is the situation where the 
Hilbert space is a tensor product Tt = H A ®Tt B of Hilbert 
spaces of two subsystems A and B. We consider finite- 
dimensional systems, and denote the dimension of TL A 
(H B ) with d-A (ds), respectively, such that the dimension 
of the tensor product Hilbert space is d = dA x ds- We 
can choose a basis of the observable algebra in A as {Ak : 
k = 1, . . . , d A } and in B as {Bk : k = 1 



basis transformation Ki = Y^,iLi ^ijMj with some ma- 
trix O then "f({Kk}) is given by 



consider the set of d\ 



d 2 B observables 



dg\, and 



{M k } = {A k ® 1, 1 <g) B k }. 



(7) 



Note that this set is not tomographically complete, since 
observables like Ak <E> Bi are missing. However, this set 
can be employed to define a very useful form of CMs. 

Definition II. 2 (Block covariance matrices). Let p be a 

state of a bi-partite system, and let M k = {Ak §5 1, 1 <8> 
Bk} be a set of observables as outlined above. Then, the 
block covariance matrix ^(p, {Mk}) has the entries 7 £j = 
(MiMj) — (Mi)(Mj) and consequently a block structure: 



A C 
C T B 



(8) 



where A — ~f(pA, {Ak}) and B = j(pB,{Bk}) are CMs 
of the reduced states of systems A and B , and 



1,3 



(AiteBj) - (A,)(B j ). 



(9) 



Similarly, we can define a symmetric block covariance 
matrix 7 S '({M / t}), for which A and B are the correspond- 
ing symmetric CMs, while C remains unchanged. 



C. Covariance matrices as description of quantum 
states 

Is it possible to completely reconstruct the state from 
a given CM? As our separability criterion uses the CM to 
decide separability, this question is important in order to 
understand, whether all states can be detected. We will 
discuss it in this subsection. Let us first show how CMs 
depend on the set of observables {Mk : k = 1, . . . , N}: 

Proposition II. 3 (Transformation of covariance matri- 
ces). Let j({M k }) be a CM as defined in (gjj. Lf {K k } 
is another set of observables, connected to the {Mk} by a 



j({K k }) = 7 ({M fe })0 2 



(10) 



Note that O is an orthogonal matrix if Ki and M t are 
orthonormal bases. 

Proof: A direct calculation gives 

(11) 



\ " 

( 



^O i<n {{M k }) l<m Ol d , 



which proves the claim. ■ 

The main point is that the previous proposition al- 
lows us to choose the basis which we want to express our 
CM in arbitrarily, since we know how the CM will be 
transformed under a basis transformation in the space of 
observables. 

We can now come back to the initial question: Suppose 
we are given some CM with a fixed basis of observables. 
Are we able to reconstruct the physical state from this 
CM uniquely? We will start answering this question by 
considering a single system. 

Proposition II. 4 (Characterization of states via 
non-symmetric covariance matrices). Given a non- 
symmetric CM with tomographically complete set of ob- 
servables, we can reconstruct the corresponding physical 
state unambiguously. 

Proof: We choose the following basis of the observ- 
ables: 



A = i = l,...,d, 



(12) 



Xij = — (|i)<j| + |j)<i|), l<i<j<d, (13) 

Y k .i = -L(|fc)(Z|_|/)(fc|), l<k<l<d. (14) 
v2 

These observables form an orthonormal basis, and we will 
refer to this basis as to the standard basis later on. As in 
any basis Mk, we can write the state as p = Ylk(Mk)Mk, 
it suffices to know the first moments (Mk). From Eq. (U|) 
one can see that 7.^ — 7^ = ([Mi, Mj]). In the following 
we will show that in the chosen basis, all first moments 
can be obtained from expectation values of commutators. 

For the chosen standard basis we can explicitly calcu- 
late all commutators 



[D k ,X, 



k,l 



V2 



Yki, [D k ,Y h 



k,l 



V2 



[X k ,i,Y k!l ] = i([k)(k\-[l)(l[). 



X k ,i, (15) 
(16) 



Hence, all expectation values of the Xi j and Y kt i can be 
calculated. The same is true for the diagonal elements: 
Using the fact that the trace of the density matrix is 
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equal to one, we can calculate all the diagonal elements 
from the mean values of [X k ,i, I 

Clearly, the same approach can be used for bipartite 
systems, if we use the CM in the full (and not in a block) 
form. In this case we can use a product basis {jti,^}}. 
Identifying (£1,12) =: i we can define the standard basis 
as above and find all first moments from the covariancc 
matrix. 

As we have seen, the non-symmetric CM defined in 
Eq. (jU describes the physical state completely. The 
knowledge of the symmetric CM in Eq. ([3]) is, however, 
not enough: 

Proposition II. 5 (Inequivalence of states and symmet- 
ric covariance matrices). The knowledge of the symmet- 
ric CM 7 s does, in general, not determine the state p 
completely. 

Proof: We prove the claim by providing a counterex- 
ample. Let us take a single qubit. As observables we take 
the appropriate normalized Pauli matrices. The symmet- 
ric CM has the following entries 



7 S - 







(loj) + fal) (!)(*,) 
4 2 



(17) 
(18) 



From this we can determine the norm of the mean value 
of the spin component in a certain direction, but not its 
sign. Hence we know the length of the Bloch vector of 
the system, up to some reflection to the origin, which 
corresponds to simultaneous change of signs of all (<Xj)'s. 

One might think that the case of one qubit constitutes 
a special case. However, the same ambiguity will arise if 
one embedded a qubit in a higher dimensional, say, three 
level system. As it can be checked, the additional observ- 
ables in the basis of observables {M k } will not provide 
any further information. ■ 

To summarize: The knowledge of the symmetric CM 
of a qubit alone is not sufficient to decide between two 
alternatives of states which have opposite (symmetric to 
the origin) Bloch vectors. Also, merely the additional 
knowledge of a single bit (the sign) is needed to make this 
correspondence unambiguous. This, however, is specific 
to the qubit case. We will now turn to investigating the 
same question for the block CM defined in Eq. ©: 

Proposition II. 6 (Relationship between bipartite states 
and block covariance matrices). For block CMs 7 and 7 
on a bipartite system, the following statements hold: 

(i) The (non- symmetric) block CM 7 determines the 
bipartite state pab completely. 

(ii) The symmetric block 7 s does not determine pab 
completely. 

Proof: Obviously, given a non-symmetric block CM for 
the set of variables Ak <S> 1 and 1 <S> Bi we can determine 



first all (Ak) and (B{) for the reduced state pa in the 
same way as in Proposition III. 41 from the blocks A and 
B of 7. Then, knowing the block C we can fix the rest 
(A k ® Bi) as 



(A k ®Bi) = C k> i + (A k )(Bi) 



(19) 



and hence (i) is proved. 

The validity of (ii) is straightforward to see for two 
qubit states, as there will be the same lack of information 
on the mean values of observables as in Proposition III. 51 
and hence 7^ B does not provide the whole information 
about the state. ■ 

The fact that the symmetric block CM 7 s does not 
determine the state completely will later be important for 
the discussion of our entanglement criteria. Therefore, 
let us investigate this correspondence for the case of two 
qubits in some more detail. For that, let Ai and Bj be 
Pauli matrices. We may write the state in the form 



PAB 



(20) 



where A,^- = tr(paf ® af ). 

As one can see from Eq. (fl"9j) we have two possibilities 
of changing the Ajj while keeping the Cjj invariant: We 
can (i) flip the signs of both of the Bloch vectors of the 
reduced density matrices, (Aoj and A^o for i,j = 1, 2, 3), 
while keeping the left hand side of Eq. ([T9]) invariant. 
Alternatively, we can (ii) flip the sign of only one of them, 
which implies that we also have to change the left hand 
side of Eq. (fT9f. 

Concerning (i), one can directly calculate the trans- 
formed state p lnv . It turns out that the eigenvalues of 
p and p lnv are the same, suggesting that they are con- 
nected by a unitary transformation maybe in addition 
with a global transposition which transforms one state 
to the other. Actually the unitary transformation is a 
local unitary one, and one has the following transforma- 
tion: 



(p inv Y =U^p AB U, 

'0 -1\ / 1 
1 8 -10 



u 



(21) 



Since there is no physical process which corresponds to a 
transposition of a state, there are two physically different 
states p and p mv which give rise to the same covariance 
matrix and which arc connected by the simultaneous flip 
of the Bloch vectors of their subsystems. Nevertheless 
we can see from Eq. (|2"Tj) that these states have the same 
entanglement properties, because there is a local unitary 
operation in addition to a global transposition connecting 
them. These transformations do not change the outcome 
of the PPT criterion, and in fact do not change the en- 
tanglement properties of any two-qubit quantum state. 

Concerning (ii) , it also possible to flip the Bloch vector 
of only one of the subsystems in a such a way that the 
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The most interesting region is the "Different" region 
where entanglement properties of the inverted state are 
different from that of the initial one. This means that 
any separability criterion which uses only the symmetric 
CM will not detect these states, as the symmetric CM is 
compatible with a separable as well as with an entangled 
state. These states will not be detected by the CMC, 
and also not by a variety of other criteria, as we will sec 
later. 



D. Properties of covariance matrices 



FIG. 1: (Color online) Entanglement properties of p £ and p" lv 
are revealed by the PPT criterion for s = 0.45 and t = -jg. e 
and r are varied. Three regions corresponds to three different 
cases. The region "Same" corresponds to the case where p e 
and p^ nv are either both separable or both entangled. The 
region "Different" corresponds to the situation where p e is 
separable but pi nv is entangled or vice versa. The last region 
consists of states p e for which the inversion of the Bloch vector 
of one of the subsystems leads to p'^ v which is not positive 
semi-definite anymore. 



whole covariance matrix will remain unaltered. This kind 
of transformation is done by 



(22) 



resulting in a transformation of p to a different p mv . Such 
a change of the state is nontrivial and can give rise to a 
matrix p lnv with negative eigenvalues, which clearly does 
not correspond to any state. Two more cases that should 
be discussed. As one can see from a numerical search, 
there are some states p, for which p lnv is still a state and 
p and p mv are either both separable or both entangled. 
But there exist also states which alter their separability 
properties after a Bloch vector inversion. As an exam- 
ple of states where p and p lnv have different separability 
properties, consider the states of the form 



/I 



V 



o 
o 

— r 
1 



(1 




p e is a slight modification of the family of the states intro- 
duced in Ref. [27j and which are known to be detected by 
PPT but not by CCNR criterion for certain parameters. 
The inverted form p™ v of this states can be calculated 
analytically. The states p £ are known to be PPT for 
(t = 0,e = 1). Going away from e = 1 and changing 
other parameters one can find regions where p e and p™ v 
have different entanglement properties. 

As we can see from Fig. [1] there arc three different re- 
gions corresponding to the different physical situations. 



In this section we will prove several properties of CMs 
which arc important for our later discussion. This con- 
cerns mainly properties of CMs for pure states and the 
behavior of CMs under the mixing of states. We will first 
show in the subsequent proposition that a suitable choice 
of observables can dramatically simplify the form of CM 
7 for pure states. 

Proposition II. 7 (Covariance matrices of pure states). 
Let Gi be a tomographically complete set of observables 
of a d- dimensional system. If p is a pure state then 7 (as 
a d x d 2 matrix) fulfills: 

(i) The rank is given by Rank(j) = d — 1. 

(ii) The nonzero eigenvalues of 7 are egual to 1, hence 
tr( 7 ) =d-l. 

(Hi) Consequently, we have 7 2 = 7- 

Proof: Without any loss of generality we assume p = 
|1)(1| and take as observables the ones of the standard 
basis (fT4|) . Calculating directly and reordering of the 
matrix elements afterwards gives a block structure [28| 



d-l 

7 = [B k ] <D d 2- 2d+2 with B k = 

k=l 



_ ( 1/2 i/2 
i/2 1/2 



(24) 

where (D^ denotes a k x k matrix of zeros. This matrix 
has the desired properties. ■ 

From this we can directly read off the properties of the 
symmetric form of the covariance matrix: 

Corollary II. 8 (Properties of symmetric CMs for pure 
states). Let {Gi} be a tomographic complete set of ob- 
servables of a d-dimensional quantum system. If p is a 
pure state, then 7 (as d 2 x d 2 symmetric matrix) fulfills: 

(i) The rank is given by Rank(j s ) = 2(d- 1). 

(ii) The nonzero eigenvalues of j s are equal to 1/2, 
hence tr(7) = d— 1. 

We now turn to a proposition concerning the trace of 
a CM for mixed states. 
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Proposition II. 9 (Trace of CMs). Let p be a mixed 
state. Then 

tr[ 7 (p)] = d-tr(p 2 ) (25) 

which implies that d— 1/d > tr( 7 (/?)) > d— 1. T/iis /lo/rfs 
also /or 7 s . 

Proo/: By definition tr( 7 ) = £ i7M = Ei^W = 
~ (Mi) 2 )- The first summation is trivial, since 
we have E/c -^fc = dl [37| . Furthermore we can write 
P = J2 k ( M k)M k which implies that £ fc (M fc ) 2 = tr (/° 2 ), 
and further 1/d < tr(p 2 ) < f. The statement for 7 s 
follows directly from the fact that tr( 7 ) = tr(7 S ). ■ 

We can also estimate the operator norm (i.e., the max- 
imal eigenvalue) of CMs. 

Proposition 11.10 (Operator norm of CMs). For the 

CM "f(p) and its linear part q(p) the operator norm is 
bounded by 

U(p)\\ < \\p\\ andhOOH < IIHI- (26) 

The same bounds hold for symmetric CMs. 

Proof: Let us first consider g(p). We have ||fl(/o)|| = 
max\ x ){x\Q(p)\x) = (a; |fl(p)|a;o) = tr(pAA r ) with 
tr(AAt) = 1. This is clearly smaller than ||p||. For 7(0) 
this follows then from (AA*) - (A) (A^) < (AA*). ■ 

Finally, CMs also satisfy an interesting majorization 
relation. This has its root in the way how one can relate 
CMs to the rotated CMs of the pure states occurring in 
their convex decompositions in terms of pure states. 

PropositionII.il (Majorization relation for CMs). For 

any (mixed) state p, both the linear part g(p) with entries 
Qi j = (MiMj) as well as the CM j(p) satisfy 

E A ^)hE A ^)] <min(M-a 73)= (27) 
j=i J=i 

for the non-increasingly ordered eigenvalues, where S 7 — 
1 for j(p) 8 1 — if q{p) is considered. 

Proof: This is a consequence of || 7 (p)||, ||fl(p)|| < 1 as 
well as of tr[ 7 (p)] < d — ^ and tr(g(p)) < d. ■ 



E. Mixing property of covariance matrices 

Separable states are those states that can be written 
as convex combinations of product states. Therefore we 
have to understand the behavior of CMs under mixing 
of states for the derivation of separability criteria. An 
important property of covariance matrices which we refer 
to as concavity property is the following: 



Proposition 11.12 (Concavity property). Let p = 

~^2kPkPk be a convex combination of states p k , then 

t(p) > ^Pk-riPk)- (28) 

k 

Clearly, this implies the same relation for the sym- 
metrized CM 7 . 

Proof: As shown in Ref . [2^ this inequality holds for an 
arbitrary symmetric CM 7 s . Moreover, since (MiMj) p = 
^2 k Pk(MiMj) Pk for all i,j, we have for the non-linear 
part that 

- (M i ) p {Mj) p > -Y^ Pk { Mi ) Pk {M,) Pk (29) 

k 

as a matrix inequality for the matrices Xij = 
-{M i ) p {M j ) p and ),. s -E k ft(«i} ft (tfj) ft . From 
this, the above inequality follows for the non-symmetric 
CM 7. ■ 

This property will later be used to derive the separa- 
bility criterion. 

F. Transformations of observables and validity of 
covariance matrices 

Transformations as generated by a general orthogonal 
matrix O used in Proposition (|IL3j) do in general not 
preserve the positivity of the state p (see Ref. for 
discussion). Only a subgroup will correspond to unitary 
transformations on the level of states. Here, we will clar- 
ify how unitary transformations of the state are reflected 
by orthogonal transformations on the level of CMs. 

For this aim, let us consider the case that p is trans- 
formed by some unitary transformation p 1— > U^pU. 
Equivalcntly, we can transform the operator basis, de- 
noted as {Gi}, as 

d ^Hi = UdU^ = E Oi,jGj. (30) 
3 

It is then easy to see that the transformation of the CM 
is 

7 (p) » 7 (p)0 T = >y(U*pU), (31) 

We can now ask in what way O depends on U , and which 
orthogonal O € 0(d 2 ) correspond to a unitary U E U{d) 
acting in state space as described above. That is, we 
look for the group representation of U (d) in the space 
of CMs (compare also the mctaplcctic representation of 
symplectic transformations in discrete Weyl systems, see 
Ref. [lf| ) . The following theorem gives an answer to this 
question. 

Proposition 11.13 (Transformation laws for CMs). Let 

U G U(d). Then the O E 0(d 2 ) representing U as de- 
scribed above i30\) is given by 

O = T T (U T <E> {/t)r*, (32) 
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where T is a d xd 2 square matrix constructed as T a ,/3\i = 
G"'^ 3 = (|Gi) , IG2), ■ ■ ■ ), where we understand a, f3 as 
a row index, G"'' 3 as vectors and construct T a ,/3\i from 
them. 

Proof: The proof is given in the Appendix. ■ 

It is an interesting open question to see how CMs 
are transformed under general completely positive maps, 
p 1 — ► AipA\, where {Ai} are Kraus operators, directly 
expressed in terms of the Kraus operators. 

At the very beginning we have given two defini- 
tions of covariance matrices for the symmetric and non- 
symmetric case (|4l5p . We will discuss this difference also 
later in the paper. However at this stage we mention a 
single connection between these two definitions for block 
CMs: 

Proposition 11.14 (Block forms of CMs under local ba- 
sis transformations). It is not possible to achieve for the 
block CMs 7 = 7 s via local basis transformations of the 
operator basis. The only states for which this relation 
holds have the reduced states pa = tr^(p) = l/^A and 
Pb = ^a(p) = where dA.B are the dimensions of 

Pa,b- It follows that 7 = j s cannot be achieved by local 
unitary operations either. 

Proof: The proof is given in the Appendix. ■ 



III. THE COVARIANCE MATRIX CRITERION 
FOR SEPARABILITY 

In this section, we introduce the covariance matrix cri- 
terion (CMC) for separability as our main topic of this 
paper. This is derived from the concavity property for 
CMs and the fact that the block CM for product states 
is block diagonal. 

Proposition III.l (Covariance matrix criterion). Let p 

be a separable state and Ai (Bi) be orthogonal observ- 
ables on TLa (Hb ), where the dimensions of the Hilbert 
spaces are dA (ds, respectively) . Define Mi = Ai ® 1 for 
i = 1, . . . ,d\ and Mi = t®Bi fori = d\ + l, . . . ,d 2 B +d\. 
Then there exist pure states |'0fc)('0fc| f or A and |</>fc)(</>fc| 
for B and convex weights pk such that if we define 

K A = Y,kPkl{\^k){tpk\) and k b = J2 k Pkl(\<i>k}(<t>k\) the 
inequality 

7 5 (P,{A^})>^®«B^ [pr b) ^ ( K Q n B ) (33) 

holds. This means that the difference between left and 
right hand side must be positive-semidefinite. If there 
are no such ka,b then the state p must be entangled. 

Proof: First note that for this special choice of Mi, for 
any product state 

l{p A ® Pb, {Mi}) = 7 (pa, {Ai}) © 70b, {Bi}) (34) 



holds. Now, since any separable state can be written as 
P = J2k Pk I ^fe) (ipk I ® \(t>k) (4>k I , we can apply Prop. (|II.12|) 
and arrive at the conclusion. ■ 

Note that the CMC is manifestly invariant under a 
change of the observables {Ak} and {Bk}, as we know 
from Proposition [TO] [see also Eq. ([57]) ]; however, a suit- 
able choice of them may simplify the evaluation a lot. 
Also, note that we have formulated the CMC for symmet- 
ric CMs, we will later discuss the case of non-symmetric 
CMs. ' 

Obviously, as such, as formulated as in Prop. IIII.ll it 
is not clear that the CMC leads to an efficient and phys- 
ically plausible test for separability: The main problem 
is to characterize the possible ka and Kb- As such, the 
formulation still contains an optimization over all pure 
product states. We will refer to an "evaluation of the 
CMC" hence whenever we can identify a property of ka,b 
that will render the above criterion an efficient test. 

Some properties of we have derived above, notably 

tr( KA ) = d A -l (35) 

(see Proposition III.8|) , which we will use subsequently. 
We will now turn to feasible ways to evaluate the CMC. 
As a first step, we have to derive conditions on the blocks 
of a block matrix as in Eq. (|33|) . which follow from the 
positivity condition in Eq. (f33|) . Then, we ask how the 
observables {Ak} and {Bk} must be chosen in order to 
make a violation of Eq. (|33|) manifest. 

Note the formal similarity of the condition 7 > ka®kb 
to tests for separability for Gaussian states for systems 
with canonical coordinates. A CM in that context |30l . l3~lj 
is any symmetric matrix satisfying 7 + ia > 0, where 

k=l V 7 

A Gaussian state p with CM 7 is now separable if and 
only if [13] there exist covariance matrices 7a, 7b, each 
satisfying 7a, 7b + io such that 

7>7a©7b- (37) 

For non-Gaussian states, a violation of this condition is 
still sufficient to detect entanglement. Eq. (|37j) . in turn, 
is simply a semi-definite program (SDP), so it can be 
efficiently decided [33[ . The characterization of the right 
hand side is here an easy task, as these matrices are again 
constrained by a semi-definite constraint. We will see in 
section lyll that for two qubits the CMC can be solved in 
a similar way as in the continuous variable case (34|. 

IV. EVALUATION OF THE CMC 

We would like to follow two strategies to evaluate the 
CMC as presented in the previous section. Both strate- 
gies are based on matrix invariants such as eigenvalues 
or singular values. Let us first characterize positive semi- 
definitcness of a block matrix of the type as in (I33|) in 
terms of singular values of its submatrices. 
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A. Evaluation of the CMC via singular values of 
submatrices 

As a start, we state the following lemma: 

Lemma IV. 1 (Block covariance matrices and unitarily 
invariant norms). If a positive matrix partitioned in block 
form is positive semi-definite, 



A C 
C T B 



> 0, 



then 



PI PI > III \c\ IP 



(38) 



(39) 



holds, where |||.||| is any unitarily invariant norm. Specifi- 
cally, this holds true for any Ky-Fan norm \\.\\kf defined 
as the sum of the largest k singular values. If we sum over 
all singular values, we arrive at the largest Ky-Fan norm, 
which is the trace norm. 

Proof: The proof of this statement is actually a corol- 
lary of Theorem 3.5.15 of Ref. [35| . It is shown that 



PI lll^l > l\(&cy/ 2 \\ 



(40) 



for any p > and any unitarily invariant norm. For p = 1 
this is the result we are interested in. We will nevertheless 
present an alternative proof of this statement for Ky-Fan 
norms ||j4||_R-F||-B||ifF > ||C|||- F , in the Appendix, as the 
proof of Proposition IIV. 9l will make use of this proof. ■ 

Using the last Lemma we have: 

Proposition IV. 2 (CMC evaluated using singular val- 
ues). Let 



A C 
C T B 



(41) 



be a CM. Then, if p is separable, we have 

\\Cf Tr < [l-tr^l)] [l-tr(p|)]. (42) 
// this inequality is violated, then p must be entangled. 

Proof: We prove the claim applying the formula (|39|) 
directly, yielding 



\\Cf Tr < \\A-K A \\ Tr \\B-K B \\ T r 



(43) 



Since A — k a as well as B — kb are Hcrmitian pos- 
itive semi-definite matrices (due to concavity property 
of CMs) their trace norm will coincide with their trace. 
Hence \\A - k a \\tt = tr(A) - tx{n A ) = 1 - tr(p 2 A ), 
where we have used Corollary 111.91 and the fact that 
E,-A< = E<((A ? ) " (^> 2 ) = (d A l) - tr(p 2 A ), since 
ti(AiAj) = Sij and p\ = J2ij(Ai)(Aj)AiAj. ■ 

Interestingly, this criterion has been proven already in 
a different context: 



Remark IV.3 (CMC and the criterion of R ef. [H[). 

The separability criterion in Proposition \IV.2\ is noth- 
ing but the separability criterion proposed in Theorem 1 
of Ref. H3/, hence the criterion of Ref. Hjj/ is a corollary 
of the CMC. 

Let us now connect the CMC to another type of entan- 
glement criteria: There are several separability criteria 
in the literature which are based on the Bloch represen- 
tation of density matrices. This representation in our 
case is just some particular choice of observables, namely 
one has to detach the identity from all others genera- 
tors, which then have to be traceless. The fact that one 
of the observables is the identity, can simplify the CMC 
sometimes. 

By definition the entries of the matrix C are given by 



C id = {A i ®B j )-{A i ){B j ), 



(44) 



which consists of a linear (in the sense of mean values) 
and quadratic part. We define £ as the linear part of C, 
i.e., £i j — (Ai ® Bj). Let us further consider £ red as the 
submatrix of £, where the entries (l A ®Bj) and (A^ 1b) 
are omitted, i.e., the first row and the first column are 
removed. Similarly, we can define matrices like 21 and 
sgrcd f rom ^4 anc [ b j n t ne same spirit, we can define 
a submatrix of k as K rcd . Note that tr(«) = tr(/t rod ), as 
the missing diagonal entry is the variance of 1, which is 
vanishing. 

We are now able to establish a connection between the 
CMC and criteria based on the Bloch representation of 
density matrices: 

Proposition IV. 4 (Relationship between CMC and cri- 
teria based on Bloch representations). Let 
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A C 
C T B 



be a CM. Then if p is separable, we have 



1 



1 



||^||| r <(!--)(!--) 



(45) 



(46) 



d A d B 

If this inequality is violated, then p must be entangled. 

Proof: First, we can define two vectors \ip) A / B with 
entries 



\iP A ) i = (A i ) |V B )i = (Si) 



(47) 



resulting in C = <£ — \ip )(ij) B \. Similar relations hold for 
A and B, so we can write the condition in CMC (f33|) in 
the form 



l* A )\ f(i> A \ 



Positivity of the left hand side implies positivity of the 
first term X alone, since we subtract only one projector 
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which is itself positive. Concerning positivity of X we can 
take » red , <B red , £ rod , and K rcd instead, since positivity of 
a matrix implies positivity of all its main minors. Using 
Eq. ([H]), we get 

||C r ° d ||Tr < P rcd - ^ d || Tr ||» rCd - K^ d \\ Tr . (49) 

Using that tr(2t rcd ) = E,> 2 <4 2 > = {d A ^A) - (l A /d A ) 
and tr(K^ d ) = (La — 1 proves the claim. ■ 

Interestingly, this separability criterion has also been 
proven before: 

Remark IV.5 (CMC and the criteri on of Ref. [16]). The 

separability criterion in Proposition \IV-4\ * s nothing but 
the separability criterion for Block representations pro- 
posed in Ref. fid] , hence the criterion of Ref. fldl l is a 
corollary of the CMC. 

Note that in Ref. [l(| the observables have been nor- 
malized in a different way, leading to a slightly different 
formula. 

Remark IV. 6 (Connection between Propositions IIV.4I 
and IIV.2[> . Proposition \IV.B is strictly stronger than 
Proposition \IV-4\ 

This fact was proven in version 5 of £13/, which came 
out after submission of our paper to the arXiv. 



and s = 1, . . . , d — 1, to conclude that 

fc 

\\A - ka\\kfm = E ^j(A - k a) 
fc fc 

= {d 2 -d+l + s)\\A\\ 

d 2 

E (51) 

j=d 2 -k+i 

where \\A\\ denotes the spectral norm of A. Using that 
ka will be more mixed in the sense of majorization than 
diag(l, . . . , 1, 0, . . . , 0) of rank d— 1 and Proposition [TT791 
one arrives at 

\\A-k a \\kfW < {d 2 -d+l + s)\\A\\-s, (52) 

and a corresponding statement for kb- Using Proposi- 
tion liV.ll one hence arrives at the observation that any 
separable state p on a bipartite Hilbert space satisfies 

((d 2 -d+l + s)\\A\\ -s) ((d 2 ~d+l + s)\\B\\ -s) 
- \\C\\ 2 KFW > 0. (53) 

It is an interesting open question whether more sophisti- 
cated uses of the knowledge of spectral properties of ka 
and Kb can be employed the further sharpen the evalua- 
tion of the CMC. 



Let us finish this subsection with a remark on the pos- 
sible use of other Ky-Fan norms in the above argument. 
In fact, we do know more about the singular values (here 
eigenvalues) of ka and kb than their sum: 

Lemma IV. 7 (Ky-Fan norms of matrices in the CMC). 
The matrices ka ( and similarly kb ) in Proposition \III.1\ 
satisfy 

fc 

E A jK0 <min(fc,d A -l), (50) 
j'=i 

for the non-increasingly ordered eigenvalues of ka ( and 
k b ). 

Proof: One can argue as in Proposition !!!. 1 1| using the 
fact that a convex combination of matrices leads to more 
mixed matrices in the sense of majorization [35| . I 

This property can immediately be applied to evaluate 
the CMC, making use of proposition IIV.1I and Weyl's 
inequalities [36[. For example, if we consider = ds = 
d and the fc-Ky- Fan norm || .||xf(*) for fc = {d 2 — d+l + s), 
we can apply the first of Weyl's inequalities with i = 1 , 



B. Evaluation of the CMC via traces of 
submatrices 

Let us first prove a simple condition on the traces of 
A, B and C, which follows from the CMC. In the fol- 
lowing, we always assume that g?a < ds- Sometimes we 
assume that the dimensions are the same, meaning that 
d = dA = ds ■ 

Proposition IV. 8 (CMC evaluated from traces). Let 

7= (54) 

be the symmetric CM of a state p and let J = 
{jii ■ ■ ■ ijd 2 } C {1, . . . ,dg} be a subset of d 2 A pairwise 
different indices. Then if p is separable, we have 

d A d A 

i=i jeJ i=i 

+ ( E S M -d B + l) = [l- tr(^)] + [1 - tr(p|)] , 
i=l 

If this inequality is violated, then p must be entangled. 
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Proof: First, note that a necessary condition for a 2 x 2 
matrix 

X = (c &) ( 56 ) 

to be positive semidefinite is that 2|c| < a + b. If p is 
separable, then by the CMC we have Y — r )~K a @i^B > 0. 
This implies that all 2 x 2 minor submatrices of Y have 
to be positive semidefinite as well. Hence for all i,j we 
have 

2|Cy < M,i + - (K A )i,i - (kb)jj- (57) 

Summing over i,j and using Corollary III. 91 proves the 
claim. ■ 

We will use this Proposition mainly for the case that 
d-A = d-B and where we sum over the diagonal entries of C. 
In this case, it just gives the condition that for separable 
states, 2tr(C) < 2 — tr(p^) — ti(p B ). This is a quadratic 
polynomial in the entries of the state, and may be viewed 
as a suitable entanglement witness on two specimens on 
a state. In the light of this fact, the criterion evaluated 
in this fashion is surprisingly strong. It is also worth 
mentioning here that one can improve Proposition IIV.8I 
by taking 4x4 minor submatrices for evaluation. Then, 
however, also off diagonal terms of ka/b will occur, for 
which not many properties are known. This makes the 
resulting conditions difficult to evaluate. 

Physically, Proposition IIV.8I says that if the correla- 
tions Ci j are sufficiently large, then p must be entan- 
gled. The question arises, how to find the observables, 
for which the Cij are large. There are several ways of 
doing this. A first result is the following: 

Proposition IV. 9 (Criterion in Proposition IIV.8I and 
diagonal block matrices). The criterion in Proposition 
\IV. 8\ detects most states if the observables are chosen in 
such a way that C is diagonal. For any state there exist a 
choice of observables that this can be achieved. However, 
even with this optimal choice of observables Proposition 
\IV.8\ delivers a strictly weaker separability criterion than 
Proposition \IV.2\ 

Proof: It is clear that the criterion is optimal, if the 
trace of C is maximal, which is the case if it is brought 
into the singular value diagonal form [2l|, ■ This can 
always be achieved [see Eq. (|97|) ]. The fact that Proposi- 
tion |IV]2] is stronger, was in a different language proven 
in Ref. [H]. 

Interestingly, the fact that Proposition IIV.8I is weaker 
than Proposition IIV.2I can also be seen from Eq. (|100|) 
from the alternative proof of Lemma TlV. II If C is cho- 
sen to be diagonal, then Proposition [TV\8] reduces to this 
equation with a = (3. Clearly, allowing a and (3 to be 
different, improves the criterion. ■ 

In the following, however, we will consider two different 
strategies: Firstly, we use the Schmidt decomposition in 
operator space of the density matrix (37j . This will lead 



to a natural choice of the observables {A^} and and 
will further connect the CMC to the CCNR criterion. 

Secondly, we will consider appropriate local filterings 
of the state [381 , l39l . l4Cl l4lj . These are active transforma- 
tions of the state, which, however, do not change the en- 
tanglement properties. Under this transformations, the 
state can be transformed into what is called its standard 
form in Ref. [Hj]. In this standard form, the CMC be- 
comes very strong and even necessary and sufficient for 
two qubits. 

C. Schmidt decomposition and the CMC 

We will first remind ourselves of what is called the 
Schmidt decomposition in operator space. It is the same 
construction as the ordinary Schmidt decomposition in 
the vector space now equipped with the Hilbcrt-Schmidt 
scalar product. A general density matrix of a composite 
system can be written as 

d\ d 2 B 

P = Y J Y<&M®G?, (58) 

k=l 1=1 

with real £k,i and the {Gf} (respectively, {Gf }) form 
an orthonormal Hermitian basis of observables. The 
Schmidt decomposition can now be achieved by diago- 
nalizing the above expression using the singular value 
decomposition of the matrix £, 

d\ 

p = Y J ^Gi®Gl, (59) 
fc=i 

where we made the assumption that dA < ds- Clearly, 
the Schmidt coefficients Afc are real and non-negative. 
Using the new basis observables {Gf} and {Gf} as ob- 
servables for the construction of the symmetric block CM, 
we have a normal form of the CMC, which we will call 
the Schmidt CMC. 

Proposition IV. 10 (Schmidt CMC). If p is separable, 
then 

2 £ I A, - >H9i9i\ < 2 - £ A [(gf) 2 + (gf) 2 ] , (60) 

i i 

where we defined gf = tr(G^) and gf = tr(Gf). If this 
condition is violated, the state must be entangled. 

Proof: Using the orthonormality of the {Gf^ B }, it is 
not difficult to see that with the observables from the 
Schmidt decomposition Cij = XiSij — XiXjg^gf holds. 
In addition, we have tr(p 2 A ) = £\ Af (gf ) 2 . Together with 
Proposition IIV. 8l this proves the claim. ■ 

Interestingly, this Proposition includes the CCNR cri- 
terion as a corollary. This shows that the CMC, even 
without filtering, and evaluated merely via the trace of 
the blocks, once the matrix is brought to Schmidt form, is 
stronger than the CCNR, which it implies as a corollary. 
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Corollary IV. 11 (CMC and CCNR). // a state p is 
separable, then in the Schmidt decomposition 



E A * 



< i 



(61) 



has to hold. This condition is just the CCNR criterion, 
hence the CCNR criterion is a corollary of the CMC. 

Proof: Using the relations \a — b\ > \a\ — \b\ and 
a 2 + b 2 > 2\ab\ we have 2 \K ~ AfoV | > 2 £, A< - 

2 Y,i x i\gfg?\ ^d 2 - E,Af[(^) 2 + ( 5 f) 2 ] < 2(1 - 

£\ \?\gfg? |), which, due to the Proposition |TV10J 
proves the claim. ■ 



D. Filtering and the CMC 

Let us now consider local filtering operations or 
SLOCC (stochastic local operations assisted by classical 
communication) [ill ] of the form 



p h-> p = (F A ® F B )p(F A ® F B 



it 



(62) 



where and F A G SL(d y i,C) and F B G SL(d B ,<D) are 
invertible matrices on the respective Hilbert spaces. 
Clearly, such operations cannot map a separable state 
into an entangled one (although they might increase en- 
tanglement measures). Also, since Fa and F B are invert- 
ible, they will also not destroy any entanglement that 
may be present in the state. In other words, p is entan- 
gled if and only if p is entangled. 

As has been shown in Refs. [H, [4l[ we can bring any 
state of full rank (i.e., p > 0) by such filtering operations 
in its standard form which is given by 



d A di 



Gi 



fc=i 



(63) 



where the {G^}, {G^} are traceless orthogonal observ- 
ablcs. Here, we again assumed that d A < d B . 

The idea now is to first apply a filtering operation and 
bring the state into its normal form. Then, the new sep- 
arability criteria are applied afterwards. Note that the 
reduction to the normal form is always possible. The 
merits of this approach are twofold: Firstly, the normal 
form reduces the number of relevant parameters, while 
still encoding all information about entanglement and 
separability. Secondly, the normal form is in a certain 
sense "more entangled" than the original state, as it was 
shown in Rcf. [io| : 

Remark IV. 12 (Extrcmality of states in normal form). 

The local filtering operations bringing a mixed state into 
its normal form are those operations which maximize all 
entanglement monotones that remain invariant under de- 
terminant 1 SLOCC operations. 



Therefore, it may be expected that many separability 
criteria become stronger if we first bring the state into its 
normal form. Note, however, that this does not hold for 
the PPT criterion, as local filtering leaves this criterion 
invariant. 

Following Ref. [3Sj . let us explain briefly an algorithm 
for transforming a state of a form as in Eq. (f55|) to its 
normal form in Eq. (|63p . As a starting point, one consid- 
ers the compact space D A <g)D B of all normalized product 
density matrices p A ® p B . For any given density matrix 
p one can define a function / of p A and p B via 



fp(pA,PB) 



tr [p{p A ® p B )] 



(detp A ) 1/dA (detpu) 



i/di 



(64) 



f P (pA,PB) is a family of positive well defined functions 
on the interior of D A ® D B , where the reduced density 
matrices both have full rank. Since p has also full rank, 
we have tr [p(p A <£> p B )\ > and because of compactness 
of D A ®D B one has even stronger tr [p(p A ® p B )\ > c p > 
0. Divergence of f p (p A , p B ) on the boundary implies that 
it has a positive minimum on the interior of D A ® D B . 

Minimization of the function f p will, as proven in Ref. 
[H[, yield the filtering operations needed. Suppose the 
minimum value for f p attained for some product density 
matrix t a ® t b with det t a > 0, dctr B > 0. Each of 
them can be decomposed as (see Eq. (66) in Ref. [H| ) 



T A T A , 



t b =T b T b , T A/B eSL(d A/B ,C) (65) 



where the T A and T B are desired local filtering opera- 
tions. Normalization factors have been ignored. 

Using this filtering operations one obtains the new 
state p which has a form 



1 



d A —1 tZ R — 1 



dA&B 



iikG\ 



G? 



(66) 



The final step involves a standard singular value decom- 
position of £ik, which leads to Eq. (|63[) . A priori, it is not 
clear whether the normal form is in some sense unique 
or not. However, it is easy to see that if we start from 
a given state and convert it into two different states in 
a normal form, then these two normal forms have to be 
connected by a local filtering operation. Using the fact 
that the reduced states of a state in the normal form 
are maximally mixed, one can further conclude that two 
different normal forms can only differ by a local unitary 
transformation. 

In practice, the minimization of f P (p A ,p B ) in Eq. (|64|) 
can be performed by an iteration as follows: let us fix p B 
and consider only the minimization over p A . This mini- 
mization can further be split into a minimization over the 
spectrum of p A and a local unitary transformation. If the 
spectrum is fixed, the optimal unitary is constructed such 
that p A and X = tr B [p(l® p B )} are diagonal in the same 
basis where the maximal eigenvalue of X has the same 
eigenvector as the minimal eigenvalue of p A and the sec- 
ond largest eigenvalue of X has the same eigenvector as 
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the second smallest eigenvalue of pa etc. If the basis is 
fixed, and (p^) are the eigenvalues of pa {X) then a 
simple calculation using Lagrange multipliers shows that 
the optimal Xk fulfill 



(y>A)/(n A « 



i^k 



1/2 



(67) 



which can be used for an iterative determination of the 
optimal Afc. In this way, the optimization can be iterated, 
converging to a minimum. Note while it is known that 
every state can be brou ght into this normal form, the 
above procedure of Ref. [3a ] is not known to be strictly 
efficient in the physical dimension d. Yet, for "reasonable 
physical dimensions" , the method in practice converges 
quickly. Moreover, and importantly, at the end of the 
procedure, one can easily (and efficiently) check via di- 
rect inspection whether the obtained filters map the state 
onto the normal form or not. Global optimality can hence 
be easily certified. 

As one can directly calculate, for a state in the normal 
form the CM takes a really simple form, namely 



1 



diag(0, d B ,d B ,. 
dAd B V dia g(°>6)6i ■ 



diag(0,fi,£ 2 ,---) 
diag(0, gU,gU,...) 



(68) 

Using this form wc obtain the following strong separabil- 
ity criterion, which we call the filter CMC. 

Proposition IV.13 (Filter CMC). If d = d A = d B and 

p is separable, then the coefficients in the filter normal 
form fulfill 



(69) 



Proof: The claim obviously follows from Proposition 
IIV.8I and the form of the CM for the normal form of the 
state. ■ 

Interestingly, for two qubits we have: 

Remark IV. 14 (Filter CMC for two qubits). For two 

qubits, the filter CMC in Provosition \IV. HH is a necessary 
and sufficient criterion for separability. 

Proof: If a two-qubit state is of full rank, the normal 
form reads 



state with finite (or vanishing) concurrence. Such a state 
will also be detected by the CMC (or not). ■ 

Direct comparison of this result with the discussion 
in Section II and Fig. [1] (and later the result of Proposi- 
tion lV.2|) might be confusing at this point, since we know 
already that the CMC itself cannot be necessary and suf- 
ficient for two qubits. This can be resolved in the follow- 
ing way: We have already learned that filtering brings 
the state in the form which in a certain sense contains 
the maximum amount of entanglement (it maximizes all 
monotones). This indeed shows that the filter CMC is 
sometimes a real improvement of the "bare" CMC, and 
filtering is more than just an appropriate choice of the 
observables. 

Let us now consider the asymmetric case, when dA < 
d B . We can formulate for this case following statement: 

Proposition IV. 15 (Separability criterion for uneven 
local dimension). If p is separable, then the following in- 
equalities hold 



d A di 



"5 



min(0, —{d B — 1) 



(4 



and 



^ [d A d B (d A - l)(d B - 1)} 



1/2 



(71) 



(72) 



holds. If one of these inequalities is violated, the state 
must be entangled. 

Proof: Eq. (|72|) is nothing but an application of Propo- 
sition IIV.4I (or IIV.2P , it has already been derived in Ref. 
[H]. Concerning Eq. dTTJ) , we will again apply Proposi- 
tion lIV.8l but with two modifications. First, when carry- 
ing out the sum over 2|Cij| < + Bjj — KA,i,i — K B,j.j 
[see Eq. (|57p in Proposition IIV.8j we do not sum over 
all Bi i. But then, we cannot subtract all of the K B ,i,i 
anymore, since d 2 B — d\ diagonal elements of k b do not 
occur in the sum. 

As a first approach, we can drop completely the sum- 
mation over all k B jj, since they are positive anyway. 
This gives 



1 3 



fc=i 



B ' 
k I > 



(70) 



where {c^ B } are the Pauli matrices (38|. Such states 
are diagonal in the Bell basis, and it is known that for 
these states Y \=i, < 2 is necessary and sufficient for 
separability [38[, 46] . Note also that the filter normal form 
can be explicitly stated for two-qubit systems. 

If an entangled (or separable) state is not of full rank, 
it can, as explicitly shown in Ref. [3^|, be brought by 
filtering operations arbitrarily close to a Bell diagonal 



d A di 



■£6 



d A 



d\ 



1 



(73) 



justifying one part of Eq. (ITTj) . 

d 2 

In a second approach, we estimate X)j=i K B.i,i- As one 
can see by direct inspection, the non-vanishing elements 
of 7 in Eq. (|68[) origin only from the linear part of CM (in 
the spirit of Proposition IIV.4I this linear part is denoted 
by g). But as we have seen in the proof of Proposition 
III. 121 that this linear part g is just the same as the linear 
part of the direct sum of ka ffi «b (denoted by Ia © ?s) 
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for separable states, i.e. g = %a®^b, hence B — <8 = ts- 
This implies that for the diagonal elements of kb the 
relation KB,i,i < = = 1/ds holds, leading to 

^ KB,i,i = d,B — 1 — KB,i,i 

> rf B -i-(4-4)_L. (74) 

This proves the second part of Eq. (fTTj) . ■ 



Corollary V.3 (Insufhciency of LUR to detect all en- 
tangled states). There exist entangled two qubit states 
which can not be detected by a LUR, hence LURs are not 
a necessary and sufficient criterion for separability. 

Proof: In the Section Hi CI we have already constructed 
a family of states p £ which cannot be detected by the 
CMC, as their symmetric block CM is compatible with a 
separable as well as an entangled state. This proves the 
claim. ■ 



V. CONNECTION TO LOCAL UNCERTAINTY 
RELATIONS 

In this section we will further analyze the connection 
of CMC with the separability criterion based on local un- 
certainty relations (LURs) [18| . To start with, we again 
state the LUR criterion as a reminder: 

Proposition V.l (Criterion based on local uncertainty 
relations). Let be A k and B k observables in system A and 
B, respectively, for which some of the variances on single 
systems is bounded by constants U a, Ub such that 

J2 > U A and J2 s2 (B k ) > U B - (75) 

k k 

Then, we have for separable states 

^8 2 {A k ®t + t®B k )>UA + U B (76) 

k 

and violation implies the presence of entanglement. 

Physically, the LURs state that separable states in- 
herit the uncertainty relations from their reduced states, 
which is not the case for entangled states. Due to this 
observation the LURs have attracted a considerable in- 
terest, and a number of interesting properties have been 
discovered: LURs can detect bound entangled states (4^ | 
and can be used to estimate the concurrence (43^ . They 
can be extended to other formulations of the uncertainty 
principle liil.lip and they can be generalized to nonlocal 
observables [29(. Finally, they can be viewed as nonlin- 
ear entanglement witnesses, which improve the CCNR 
criterion [371 ]. 

For the connection to the CMC we have the following: 

Proposition V.2 (Connection to local uncertainty rela- 
tions). A state p violates the CMC for symmetric GMs 
iff it can be detected by a LUR. 

Proof: The proof is given in the Appendix. ■ 

This result show that the LURs for appropriate observ- 
ables and the CMC are equivalent, however, the CMC has 
the major advantage that it can be directly evaluated, 
while for the LURs the appropriate observables have to 
be identified. Moreover, we can state: 



VI. THE CMC FOR TWO QUBITS 

After the previous discussion of the situation of Hilbcrt 
spaces of arbitrary finite dimension, we now turn to the 
important simple case of a 2 x 2-system - two qubits - in 
some more detail. We take as observables the set {A k } = 
{B k } = {l/V2,a x /V2,o- y /V2,o- z /V2} as in Eq. Q. 

Since these observables contain the identity, one can 
easily check that many terms in the symmetric block CM 
vanish. Effectively, 7 is actually a 6 x 6 matrix (denoted 
by 7° ff ) originating only from the {A k } and {B k } with 
k = 1,2,3 which are not proportional to the identity, 
and not by an 8 x 8 as one could guess from the general 
theory 

To characterize the ka in the CMC, note that for a 
pure state |a) on system A we find, according to Propo- 
sition III.8I the following properties of the 4x4 matrix 
7(|o)<o|): 

(i) Rank(j) = 2. 

(ii) The nonzero eigenvalues of 7 are equal to 1/2 in a 
suitable basis. 

We also know that in the chosen basis, the first row as 
well as the first column of 7(|a)(a|) vanish, and we have 

7(|a)(a|) = © 1 e7(|a)(a|r ff (77) 

where 7 is the effective 3x3 CM as above. This has 
to be of rank two with eigenvalues 1/2. This implies that 
7(|a)(a|) can be written as 

7(l«X«ir ff = ^(l3-|0a)(0 a |), (78) 

where 1 3 denotes a 3 x 3 identity matrix, and \<j) a ) € R 3 . 
In fact, any matrix of this form is a valid CM: 

Lemma VI. 1. For any vector \<f>) € R 3 a matrix of the 
form (I3 — \4> a ){<p a \)/2, is a valid CM of some two qubit 
state. Consequently, the set of valid ka is given by all 
matrices of the form 

*a = \(U-Pa) (79) 

where pa is a real 3x3 matrix with trace one and positive 
eigenvalues. 
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Proof: We have already shown that the CMs are of the 
required form, and only have to argue that any matrix 
of the form X = (1 3 - \(j) x ){(j) x \) /2 is a valid CM. To see 
this, note that unitary transformations of the \a) result in 
orthogonal transformation on j(\a)(a\) . Moreover, for 
the special case of a single qubit any orthogonal trans- 
formation on 7 eff can be generated by a unitary transfor- 
mation on state space [46j . expressing the isomorphism 
between the Lie-algebras su{2) and so(3). Therefore, we 
can transform X into 7(|a)(a|) off and construct the cor- 
responding state vector \x). To finish the argument, note 
that the set of all ka is by definition the set of all convex 
combinations of pure state CMs. ■ 

After having proven this Lemma we can formulate the 
two-qubit version of the CMC: 

Proposition VI. 2 (CMC for two qubits). Let p be a 

state of two qubits and let 

{A k } = {B k } = {t/V2,o-yV2,ayV2,o- z /V2} (80) 

be the chosen set of observables. Let be the 6x6 CM 
as mentioned before. Then the state p fulfills the CMC 
iff there exist 3x3 density matrices pa and ps such that 



\ls + \{pA®p B )>Q- 



(81) 



Proof: The claim follows if we insert the k's from Eq. 
(|79f into Proposition IIII. ll Note that it suffices to find 
complex pa and ps- If we can identify such matrices, 
their real part will saturate Eq. (|8"Tj) as well. ■ 

In this form, the problem is a special instance of an 
efficiently solvable semidefinite program (SDP) [33[ in 
primal form, a feasibility problem . 

In general, a SDP consists of a linear function c T x 
which is minimized subject to a semi-definite constraint 



F(x)=F + J2^i F i>0, 



(82) 



which is linear in the problem variables Xi. Hence the 
problem is defined by the real vector c and by the hermi- 
tian or symmetric matrices Fi. If c — 0, then the prob- 
lem is referred to as a feasibility problem. Via Lagrangc- 
duality, a dual problem can be formulated in which the 
expression —ti(FoZ) is maximized over a positive semi- 
definite (hcrmitian or symmetric) matrix Z, with the con- 
straints that tr(FiZ) = Ci. Since 



c T x + tr(F Z) = tr(F(x)Z) > 



(83) 



holds true due to the positive semi-definiteness of F(x) 
and Z, solutions of the dual problem deliver a bound 
on the solutions of the primal problem and vice versa, 
which is referred to as weak duality. Finally, if there is 
a solution to the primal problem with F(x) > or a 
solution to the dual problem with Z > 0, then strong 



duality holds, meaning that a pair (x* , Z*) exists such 
that c T x* + tr(F Z*) = holds. See also Ref. H3 for an 
extensive treatment of the subject. 

For the evaluation of the CMC, we can formulate the 
problem differently, such that if the primal problem de- 
tects the state as entangled, then from the solution of 
the dual problem local operators can be extracted which 
allow for the detection of the state with LURs. This is 
similar in spirit as the solution in the continuous variable 
case 



34 



Explicitly, we formulate the primal problem as 



min —A 

„eff 



(84) 



subject to 7° — ka ffi «b > 
1 



k a ,b = -[(1 + A)1 3 -/9a,b] >0 

tr(pA,s) = 1 + A. 

In this formulation, the matrices ka,b ar e positive and 
have trace 1 + A. If the constraints can be fulfilled for 
A < only, then the state corresponding to 7 cff is entan- 
gled. The SDP can be formulated with block-diagonal 
matrices {Fi} collecting all the constraints. For instance, 
by inserting the definition of Ka,b into the first constraint 
and expressing the equality constraints by a '>' and a '<' 
constraint, we can write Fq as 



1 



1. 



1 



Fq = (7 -2 l6 )®2 l30 2 l3ffi( " 1)el ® ( " 1)e1 ' (85) 

and the matrices Fi accordingly by choosing a basis for 
real, symmetric matrices for the blocks. Without loss 
of generality, the matrix Z can be chosen block-diagonal 
accordingly. In the order from above we have Z = Z\ ffi 
i Z* 1 ffi Z£ 2 © Z^ 1 ffi Z? 2 , where Z x is a 6 x 6 



Z*< 



A B 

matrix, Z 2 ' are of dimension 3x3, and Z. 
scalar. The dual problem can then be formulated as 



A,B;1,2 



are 



max 



[tr( 7 et %) - 1] 



(86) 



subject to --[tv(Zx) - tr(Z^) - tr(Zf 



zf - zf + zi 1 - zi 2 - 1 

I rrJ\,a\ I rrA,B\ r./ryA.B.l yA,B;2\ 



A,B\ 
1 > 
A,B\ 



(ZD 



><J 



(Z 2 ' )i<j> 



A B 

where Z 1 ' are the single-particle subblocks of system 
A and B, respectively, and i and j run from 1 to 3. It 
turns out that Z\ has the properties of an entanglement 
witness in the space of covariance matrices (CM-witness) 
as in the continuous- variables case [34[ : 

Proposition VI. 3 (CM- Witness from dual program). 

For every feasible solution Z to the dual problem formu- 
lated above, the matrix Z\ is a CM-witness in the sense 
that it fulfills tr(jg^Zi) > 1 for all CMs 7^ from sepa- 
rable states. Hence if Z\) < 1 then the correspond- 
ing state is entangled. Further, it is optimal in the sense 
that tr(7 e ^Z 1 ) is the minimal value of tv{"f e ^X) for any 
X > of the same dimensions. 
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Proof: It follows from weak duality that tr(7 eff Zi) > 
1+A, hence tr(^g S Zi) > 1 holds for all 75 s from separable 
states. In this case, strong duality holds, which we prove 
by providing an example: 



■If 



Is 



3 
4 



1 > 



(87) 



fulfills all constraints. Hence there exist (X* , Z*) such 
that tr(7 off Z*) = 1 + A* holds, and the dual program 
reaches the minimal value of tr(7 eff Zi). ■ 

If the entanglement of a state is detected by a CW- 
witness Z\, then it is possible to write down a LUR de- 
tecting the state as well. This is remarkable because it 
is in general very difficult to find a LUR detecting the 
entanglement of a given state. 

Proposition VI. 4 (LUR observables from witness). 

Given a CM-witness Z\, it is possible to define LUR ma- 
trices {Ak} and {B k } from Z\ such that 



tr(7 e -%i) = J2 s2 (Ak ® 1 + 1 ® B k ) 



(88) 



holds. 



Proof: The LUR corresponding to Z\ can be extracted 
as shown in the proof of Proposition lV.2l in the Appendix: 
we can spectrally decompose Z\ = J2 k •^fe|'0fc)('0fe| =: 
Efc A fe |c*W © /3( fc ))(a( fe ) © /3( fe '|. Defining the local LUR 

variables A k = y/X~kJ2i ^ A i and B k = V^~kJ2i P[ k) Bi 
we have for p that tr(Zi7° ff ) = £ fc S 2 (A k © 1 + 1 © B k ), 
where {A k } and {B k } are defined in Eq. ([80]) . ■ 



VII. EXAMPLES 

In this section, we consider bound entangled states of 
two different types, and investigate the strength of the 
different criteria discussed in this paper. 

In the first example, we take the 3x3 bound entangled 
states, called chessboard states, introduced by D. Brufi 
and A. Peres (48|. They are defined as 



P = AT y £\V j )(V j \, 



(89) 



where Af denotes the normalization, and we used the 
unnormalizcd vectors 

|Ui) = \m,0,ac/n;0,n,0; 0,0,0), 

\v 2 ) = 10,0,056, 0,050,0,0), 

|V 3 ) = |n,0,0;0,-m,0;ad/m,0,0), 

\Vi) = |0,6 ) 0;-o,0,0;0,d,0). (90) 

Characterization of the family is done by six real param- 
eters. We tested all criteria, presented in this paper on 
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FIG. 2: (Color online) Detection of 3 x 3 chessboard states. 
For the different criteria the fraction of states which are de- 
tected is shown. See text for further details. 



randomly generated chessboard states, where parameters 
have been drawn from the normal distribution with zero 
mean value and standard deviation of two. The results 
of this test are presented on the Fig. [21 

As one can see from Fig. [2] the most of the states are 
detected by bringing first the state in its normal form 
(Proposition IIV.13|) - 98.86% of all states. The criterion, 
which uses an estimation of singular values of the off di- 
agonal block of CM (Proposition IIV.2[) , which was also 
proposed earlier in [2fj| detects 22.57%, whereas another 
criterion proposed in this paper (Proposition IIV.9[) de- 
tects 22.00%. Moreover the criterion, which uses Schmidt 
decomposition (Proposition [TVJjB detects 20.00% which 
is more or less the same amount as is detected by CCNR 
criterion - 19.52%. Finally the criterion presented in 
Proposition lIV.4l which was first proposed by de Vicente 
fl6| detects only 8.57% of randomly generated chessboard 
states. 

As the second example, we consider 3x3 bound en- 
tangled states arising from an unextendible product basis 
14911 . mixed with white noise: 



l^o) 
1^2) 
1^4) 

PBE 

Pup(p) 



1 



V2 



|0)(|0)-|i», bAi 



1 



(|o>- 



^=|2>(|l>-|2», 



I<|°> 



|2»(|0) + 

4 



= V! (|1> 

|1> + |2», 



|1»|2>, 
|2»|0>, 



PPbe 



(91) 



These states are detected by Proposition IIV.13I for p > 
0.8723 while the best known positive map detects them 
only for p > 0.8744 (see [37} and references therein). Be- 
sides this we have also tested all other criteria presented 
in this paper. Criteria of Propositions lIV.2l IIV.9l both de- 
tect these states for p > 0.8822. The criterion derived for 
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Schmidt decomposed states (Proposition IIV. 10| detects 
the states for p > 0.8834, whereas the CCNR criterion 
detects them for p > 0.8897. Finally Proposition IIV. 41 
detects the states for p > 0.9493. 

Finally, let us shortly comment on the efficiency of the 
implementation of all these criteria. The filtering oper- 
ation can be implemented quite fast, using the simple 
algorithm outlined above takes a few seconds on a desk- 
top computer (5x5 system: ca. 6 sec, 10 x 10 system: 
ca. 24 sec, 15 x 15 system: ca. 72 sec). Then, the trace 
norm of C can be quickly computed as the trace norm of 
the realignment of the matrix p — pa® Pb [SJ . For com- 
parison, only the first step of the semidefinite program of 
Ref. [l3[ requires already ca. 10 min. for a 4 x 4 system, 
becoming practically unfeasible for higher dimensions. 

VIII. CONCLUSION AND OUTLOOK 

In this work, we have further developed the ideas of 
Ref. [22| and investigated the covariance matrix criterion 
(CMC). We have shown that this is a strong separabil- 
ity criterion, which can be simply evaluated. Combined 
with filtering it is necessary and sufficient for two qubits 
and in higher dimensions it detects states where the PPT 
criterion fails. Moreover, it contains many other separa- 
bility criteria, which have been proposed to complement 
the PPT criterion as corollaries. 

There are several open problems which deserve a fur- 
ther investigation: 

• First, one might study the exact relation between 
the CMC for symmetric CMs and non-symmetric 
CMs. In the present paper, we have used only the 
trace of n A l B for the evaluation, hence this differ- 
ence did not become apparent. However, as the 
non-symmetric CM describes the state completely 
and encodes therefore all information about the 
separability properties, this difference might be a 
way to improve the CMC. In addition, one could 
investigate the relation between the linear part and 
the nonlinear part of the CM in some more detail. 
In the proof of Prop. IIV. 151 we have seen already 
that such an investigation may indeed improve the 
CMC. 

• Another interesting open question is to relate the 
CMC to quantitative statements, such as to the 
estimation of entanglement measures. One should 
expect that the program of Ref. [Fl| linking vio- 
lations to quantitative statements on the entangle- 
ment content is applicable to the discussed criteria. 
Also, similar relations hold for Gaussian states • 
First ste ps i n this direction have already been taken 
in Refs. [13, llllil- 

• Finally, it would be very interesting to develop a 
theory similar to ours for entanglement of multipar- 
ticle systems. Here, however, a significant amount 



of work has yet to be done, as it is not even obvi- 
ous how to identify the object corresponding to the 
block CM for multipartite systems. 
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X. APPENDIX 

In this appendix, we present the more technical proofs 
of our previous statements. 

Proof of Proposition I II. i 31 Let us first explain some 
properties of the matrix T. This matrix has entries which 
are just the basis vectors G, written as columns. More- 
over, rrt = 1 = r+r, i.e., r is a unitary, since 

fe a. .p 

= Y^^ i ' a )Gf = tv{G i G j )=5 i<j , (92) 

a. .[3 

where we have used the orthogonality and hcrmiticity of 
Gj. However Y is a special unitary, since the columns 
correspond to orthonormal Hermitian observables. Now 
we have in Eq. (|30jl 

Y. (: ';'( <)I X., c^' (93) 

where we have used the definition of F and the fact that 
the expression in the middle of Eq. (|93[) is nothing but 
i-th column of TO T . Conversely, 

= (U® U*) a ,p, s „ I* i7 n = (U® U*T) aM , (94) 

where we have used the definition of T and that A^k <& 
Bi m = (A® B)ij t k,m- Therefore we can write 

T = T ] (U ® U*)Y = T T (U* ® U)T*, 

o = r t (?/ t <S)U T )r = r T (u T <z>u^)r*, (95) 

where we used that O is real. With these representations, 
we can finally check the orthogonality of the O as 

O t O = T T {U* ®U)T*T T (U T <x>[/ f )r* 
= v T {u* ®u)t{u T <x>[/ f )r* 
= v T {u*u T ®uu v )T* = r T r* = 1. (96) 
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Proof of Proposition \II. 14\ First note that if we write 
7 as in Eq. (JSj) in the blockwise form, A, B correspond 
to CMs of the subsystems A, B and C has entries of the 
form (Ai © Bj) — (A4) (Bj), where Ai, Bj are observables 
taken for subsystems A, B. 

The condition 7 = 7 s is equivalent to the condition 
7 = 7 T , in particular A = A T and B = B T . If we change 
the local bases on A and B via O — O a © O b the CM 
gets to 



B\T 



7' = (o A © O b )~/{O a © O a ) 



(97) 



As we can immediately see -y' T = (O a © O b )^ t (O a © 
B ) T = 7' if and only if 7 T = 7, so the symmetry of 
CM does not depend on the particular choice of basis in 
observable space. 

Therefore we are able to choose the standard basis. 
Let us consider only subsystem A, i.e., left upper block 
of matrix 7, and let us assume that A = A T holds. As we 
have showed already, we can obtain pa from the matrix 
A by use of the commutators Aij - A hi = ({M A , M A }). 
However, all those commutators vanish for the case A = 
A T . Since then (Xk,i) = (¥k,i) = for all k, I, it follows 
that pa is diagonal. The diagonal elements can be also 
determined as in Prop. HL4l and since also (Z k j) = for 
all k,l, it follows that pa = 1/d.A, which completes the 
first part of the proof. 

Finally, local unitary transformations are only a sub- 
class of the orthogonal transformations considered be- 
fore, hence 7 = j T cannot be achieved by a local unitary 
transformation of p neither. ■ 



Alternative proof of Proposition \IV.1\ for Ky-Fan 
norms: For a matrix as in Eq. (|38|) the following con- 
dition has to be fulfilled: 



A 
C T 



(l«> l/3»>0, 



for all vectors \a), \f3), which implies that (a|^4|a) + 
(P\B\/3) >2{a\C\f3), where we took -\0) instead of \0) 
for convenience. Especially, we can take \a) = a\ip k ) and 
|/3) = /3|</>fc), where the vectors \ip k ) and \4> k ) are singular 
vectors from the singular value decomposition of C and 
(ipk\C\(j>k) = crfe(C) is the fc-th singular value. Hence 

a 2 {^ k \A\^j k ) +(3 2 {cj) k \B\cj) k ) > 2ap{il> k \C\<l>k). (99) 

Note that (tfj k \A\i(j k ) and (tfj k \A\i(j k ) are greater than zero, 
because A and B are positive semi-definite matrices. 
Taking the sum over k and noting that for A and B ex- 
pressions like y*£_i (il>k\A\il}k) are a lower bound on the 
K-th Ky-Fan norm [35[ we get 



KF 



KF 



> 2a0\\C\\ 



KF- 



(100) 



The last formula is necessary and sufficient condition for 
the 2x2 matrix 



\M\kf 
\C\\kf 



\C\\kf 
\B\\kf 



> 



(101) 



to be positive semi-definite and having a non-negative 
determinant, from which the claim follows. ■ 



Proof of Proposition I V. 2k The proof is an adaption of 
a similar proof given in Ref. [2i|. We will often use the 
property that CMs can be used to compute variances. 
Imagine TV = v^M^ is a linear combination of the 
M fe with v k G R, then 

$ 2 {N) = Y J Vkn{{M k M l ) - (M k ){Mi)) = (v\i({M})\v). 

k,l 

(102) 

Let us now assume that p violates the LURs and we 
can find A k , B k , Ua and Ub as in Proposition IV.ll Wo 
assume that the CMC is fulfilled, i.e., there exist ka and 
kb such that for the CM 7 we have 7 > ka © «s- We 
can write 



A k = Y, a( i k)Ai and = A (fe)fi i 



(103) 



where the {A k } and {B k } are the observables chosen in 
the definition of 7. This leads to S 2 (A k © 1 + 1 © B k ) = 
(aW © /3W|7|aW © (3^). Also, by definition 

KA © KB = 5>i7(|ai><ai|)©7(l&J><&i|) (104) 



and hence (a.W © © K B \a^ © = 

X; z P/[<5 2 (-4fe)|a i )<Q ! | + S 2 (B k )\ bi){bi\]- But then summing 
over k yields 



fc 

>Ypi [ s2 ( A k)\ ai )( ai \ +S 2 (B k )\b l){bl \] 



(105) 



k .i 



> 



\a){a 



k)\b){b\ 



> Ua + Ub, 



which is a contradiction to our assumption that p violates 
the LURs. 

To show the converse direction, let us assume that 
q violates the CMC. Let us define a set of matrices as 
X = {x\x = ka® kb + P with P > 0}, which geomet- 
rically is a closed convex cone. Using this definition, we 
can formulate the CMC differently, by saying that if p is 
separable, then 7 G X. As our p violates the CMC, we 
have 7 ^ X. 

According to a corollary to the Hahn-Banach theorem 
[50| for each 7 ^ X there exist a symmetric matrix W 
and a number C such that tr(W7) < C while 



tr(Wx) > CVx e X. 



(106) 



Since X is a non-compact cone, and we can add arbi- 
trary positive operators to the elements of X, we can 
conclude that ti(WP) > has to hold for all P > 0, 
and consequently we have W > 0. Now let us make 
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use of spectral decomposition of W and write W = 
EfcAfclVkWfcl =:£ fc ^|a (fc) ©/3 (fc) }(a (fc) ©/3 (fc) l- Defin- 



ing A k = y0^kJ2l a l A l aild B k = V^fcE/A B l WC 

have for p that 

tr(U/ 7 ) = ^5 2 (i fe © l + l® B fc ). (107) 

Furthermore, by definition we have for all k a © kb € X 
and from Proposition III. 121 it follows that all 7 a © 7b G 
X. Hence for each product state p = pa © Pb we have 



C < tv{W lA © 7s) 
implies that 



E k [5 2 (A k ) PA +SHB k ) PB ]. This 



C< mm [V(<5 2 (i,) P4 +<5 2 (S 



k) pi 



Finally, since the CMC is violated, 7 ^ 
J2 k <5 2 (i fe © 1 + 1 © B fc ) - trQVi) < C < 
leading to a violation of the LURs criterion. 



X 
U A H 



and 

U B 



Note that in principle this proof also applies to the 
CMC for non-symmetric CMs. Then, however, the "ob- 
servables" in the LURs will be non-hermitian, their vari- 
ance has to be defined as S 2 (X) = (XX^) - (X){X^) and 
their physical interpretation is not so clear. ■ 



<min[J2S 2 (A k ) PA ] +min[J2S 2 (B k ) PB ] 

PA PB 



U a + Ub 



(108) 
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